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Abstract. Two strongly coupled quantum dots are theoretically and experimentally 
investigated. In the conductance measurements of a GaAs based low-dimensional system 
additional features to the Coulomb blockade have been detected at low temperatures. These 
regions of finite conductivity are compared with theoretical investigations of a strongly 
coupled quantum dot system and good agreement of the theoretical and the experimental 
results has been found. 
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1. Introduction 

The electron transport through a single quantum dot in a two-terminal configuration is 
governed by the interaction energy of the electrons on the dot. In many experiments, the 
state of the quantum dot is essentially characterised by the electron number since orbital 
degrees of freedom do not play a major role and can thus be ignored. In the limit of weak 
dot-lead coupling, the resulting current is determined by states with an energy above the 
Fermi energy of one, but below the Fermi energy of the other lead. The other states suffer 
Coulomb blockade: Energy conservation together with Pauli's exclusion principle preserves 
their occupation number and, consequently, they cannot contribute to the transport. Only 
when the dot- lead coupling becomes larger, co-tunneling processes start to play a role and 
suspend Coulomb blockade. 

When two or more quantum dots are in a linear transport arrangement between two leads, 
the inter-dot tunneling can be incoherent or coherent, depending on the coupling strength. 
Incoherent tunneling is sequential, i.e. between two tunneling events, the electrons dwell in 
one particular dot. Coherent tunneling is found for strong inter-dot coupling such that the 
electrons reside in the delocalised eigenstates of the double dot. The analogy to ^-electrons 
in molecules is reflected by the term "artificial molecule". A convenient theoretical picture 
for coherently coupled quantum dots is a one single central system in which orbital degrees 
of freedom play a role. 

In an unbiased double dot, the relevant orbitals are the bonding and the anti-bonding 
superposition of the localised states. Then an electron prepared in one dot will tunnel forth and 
back to the other dot with a frequency set by the tunnel splitting. These coherent oscillations 
can be observed by lowering after a waiting time the chemical potential of, say, the right lead. 
If at that stage the electron is in the right dot, it will tunnel to the right lead. Periodic repetition 
of this procedure yields a dc current that reflects the coherent oscillations [1]. The coherence 
of the superposition together with the possibility to perform a readout allows one to devise 
charge qubits with double quantum dots. The orbital degrees also influence the transport under 
microwave excitation: Microwave irradiation can induce electron transitions from the ground 
state to an excited state and thereby enhance the electron transport between the leads, so that 
one observes photon-assisted tunneling [2-4]. 

A further common method to characterise low-dimensional semiconductor systems such 
as the mentioned quantum dots, are conductance measurements at low temperature: Since the 
current changes whenever an energy level enters or leaves the voltage window, the differential 
conductance exhibits a corresponding peak. Shifting, in addition, the energy levels by a gate 
voltage yields the characteristic "Coulomb diamonds" which are observed in the differential 
conductance as a function of gate voltage and bias voltage. Within this work, we study 
both theoretically and experimentally the fingerprints of orbital degrees of freedom in the 
Coulomb diamond structure of coherently coupled quantum dots. In section 2, we describe 
our experimental setup and present transport measurements, while in section 3, we study 
a minimal model that exhibits the observed Coulomb diamonds. Moreover, we relate our 
theoretical findings to the experimental data. 
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Figure 1. (a) Scanning electron micrograph of the sample. Red (darker) areas are highlighting 
2DEG, black (brighter) areas are wet-etched and non-conducting. The annealed ohmic 
contacts (not shown) are labelled. The depletion of the left channel (under investigation) is 
controlled by the left in-plane gate. The channel can be made narrower by applying a more 
negative voltage to the gate until all electrons are forced out of the channel. In our experiment 
the right channel is completely depleted by a applying a sufficiently large voltage (— 4.5 V) to 
the right gate as shown in panel (b). (c) Wiring scheme of the sample. 



2. Experimental setup 

In this work, we used a GaAs/AlGaAs heterostructure for sample fabrication, where the 
two dimensional electron gas (2DEG) is located approx. 55 nm beneath the sample surface 
[5]. First, a mesa structure was defined by photolithography, followed by wet-etching and 
annealing of ohmic contacts for source, drain, and in-plane gates. Then the nanostructure 
formation has been processed by electron beam lithography and wet-etching [6]. Figure 1 
shows the nanostructure under investigation. While the grey (light) areas were etched, the 
red (dark) areas depict the regions containing a high mobility 2DEG. The structure is 8 /im 
long and in total about 10/im wide. The inner structure is sawtooth shaped with four teeth 
and asymmetric with respect to the vertical centre line. The conductive channels between the 
sawtooth-tip and in-plane gates are approx. 0.9 /im wide. The 2DEG has an electron density 
of 3.95 x 10 15 /l m 2 an d a mobility of 51.7m 2 /Vs (both measured in the dark at T = 4.2K). 
The sample allows us to carry out measurements individually on each channel. Here, we 
report only on measurements on the left channel. In order to ensure that only the left channel 
is conductive and the measurement is not affected by the right channel, a relatively high 
negative voltage of —4.5 V is applied to the right in-plane gate. This causes a depletion [7, 8] 
of the 2DEG in the right channel such that no electrons can pass from source to drain. This 
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Figure 2. Measured current through the left channel with decreasing voltage on the left in- 
plane gate at 1 .3 K. The applied source-drain- voltage is Vsd = 0.1 mV. The system is evolving 
from the Id regime, characterised by current plateaus, the into QD-regime with characteristic 
Coulomb oscillations. Inset: blow-up of the very first oscillation ([/gate = —2.66 V) next to the 
pinch-off. 



situation is shown in figure 1(b), while figure 1(c) sketches the wiring scheme of the sample. 
In the same manner, the left channel can be depleted as well. If a negative gate voltage is 
applied, the channel becomes narrower until three potential barriers between the depleted 
areas, developed from the voltage on the gate and the wet-etched tines, are formed. Now, the 
Id channel is separated into shorter channels and with decreasing gate voltage it eventually 
evolves into small quantum dots. At a certain voltage, the so-called pinch-off, the channel 
is completely depleted and, consequently, no current can flow. Such an evolution from ld- 
channel [9, 10] over quantum dots (QD) [10-14] to total depletion is shown in figure 2. At less 
negative gate voltages, characteristic ld-conductance quantisation in the form of current steps 
— conduction plateaus — can be seen. At even more negative voltages, current oscillations 
were observed, being characteristic for QD transport. The irregularity of the current peaks 
as a function of the side gate voltage already indicates the existence of a rather complicated 
electronic structure close to the pinch-off. The measurement in figure 2 was carried out at 
temperature T = 1.3 K with source-drain bias voltage of VsdO.1 mV. 

Because it is quite difficult to estimate the exact landscape of the barriers nor to explain 
the measurements in greater detail if the channel consists of more than one quantum dot, 
we focused on measurements very close to the pinch-off, where the first current oscillation 
appears. The first current peak is at Vg a t e = —2.66V; see inset of figure 2. This quantum 
dot system can be characterised by a set of current measurements for different source- 
drain voltages [15]. Figure 3 shows the corresponding differential conductance d//dVsD 
which exhibiting characteristic Coulomb-diamond structure [11-13, 15]. From the slopes 
of the transition from high to low conductivity shown in figure 4, one can extract the 
effective parameters which we use later in our theoretical description [15]: the capacities 
C g ate = (4± 1) x 10~ 18 F, C drain = (51 ± 15) x 10~ 18 F and C z = (103 ±8) x 10~ 18 F where 
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Cgate is the capacity between the dot and the gate, Cd ra in the capacity dot-drain, and Cz the 
total dot capacity. The single-electron charging energy is Eq = 0.78±0.06meV and the 
energy spacing between two levels inside the dot is AE = 1.6 ±0.1 meV. The relatively large 
uncertainties follow from the estimated uncertainty of reading the quantities from figure 3 
and from the calculated propagation of uncertainty. If one assumes a parabolic potential, the 
energy spacing lets us estimate the lateral dimension of (14± 1) x 10 _15 m 2 of the dot and the 
diameter of ( 136 ± 6) nm. However, because of the triangular shape from the lithographically 
defined sawtooth potential, it is unlikely that the dot is perfectly round and a irregular shape 
is assumed. 

Figure 3(a) shows the differential conductance d//dVsD- It exhibits some differences as 
compared to "regular" Coulomb diamonds. First, the diamond is slightly tilted to the left. This 
is an indication for asymmetric tunnelling barriers between the quantum dot and the leads as 
seen in almost every Coulomb-diamond measurement. Second, some additional structures in 
the diamond can be spotted. The most obvious one is the small area of high conductivity in the 
center section where the tips of conductive areas almost merge (marked with * in figure 4). A 
further interesting feature is the narrow stripes of finite conductance alongside the main areas 
(marked with +). These mentioned areas are also rather symmetric due to bias source-drain 
voltage and asymmetric along the gate voltage, i.e. they appear only for lower gate voltage 
whereas the transition from electron transport to Coulomb blockade for higher gate voltages is 
very sharp. The third area (marked with #) exhibits a high conductivity and is also symmetric 
in the source-drain voltage but not as a function of the gate voltage. 

3. Theoretical descripton 

For a master equation description of electron transport at very low bias voltages, one needs 
to take particular care in order to avoid inconsistencies like the emergence of spurious 
non-vanishing transport in equilibrium situations. Such problems typically arise from the 
approximation in the interaction representation of the coupling operator [4, 16]. An detailed 
derivation of such a master equation approach has been presented e.g. in reference [17]. Here, 
we will briefly review this approach. 

3.1. Model 

The setup at hand for studying coherent quantum transport is shown in figure 5 and the 
corresponding Hamiltonian reads 

H = //dots + Pleads + //dots-lead • ( 1 ) 

The individual terms describe the quantum dots, the electron reservoirs of the leads, and the 
coupling of the dots to the leads. The system itself is treated in a tight-binding approximation 
which we restrict in the following to two orbital degrees of freedom. Since we aim at exploring 
blocking effects, the corresponding wire Hamiltonian, incorporating the Coulomb repulsion 
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Figure 3. False-colour-plot of the differential conductance d//dVso as a function of source- 
drain bias and gate voltage. Blue corresponds to low conductance and red to high conductance. 
The experimental data (a) described in section 2 are compared to theoretical results for a double 
quantum dot with interacting (b) and non-interacting (c) electrons. The theoretical calculations 
are for dot-lead couplings Tl = 0.2A, Tr = 0.25A and temperature T = 0.1 A. 



in the limit of a large interaction strength U, assumes the form 

//dots = ££„cJc„-A (cla +c\c 2 ) +UJS(J/ _ i). ( 2 ) 

n 

The fermion operators c\ (c n ) create (annihilate) an electron in the orbital |n) and E n denotes 
the respective on-site energy. In the Coulomb interaction term, ,jV = c\c n is the operator 
counting the excess electrons on the dots. The inter-dot coupling is characterised by the 
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Figure 4. Sketch of the system characteristics superimposed to the experimental data.From the 
slope of the dashed lines one can estimate the capacities of the dot, the doted lines highlight the 
stripes of elevated conductivity. The symbols *, + and # mark parameter regions mentioned in 
the text. The total capacity of the system, Ci, can be estimated from the value of AVsd [15]. 




Figure 5. Tight-binding model for a double quantum dot coupled to two leads. An external 
bias voltage V = (/Ir — jUl)/<? is applied to the mesoscopic system. 

hopping matrix element A. The leads attached to the dots are modelled by ideal Fermi gases, 

#leads = J! Yl e Q C £q Ce V ^ 

£=L,R q 

where c\ q (c£ q ) creates (annihilates) an electron with energy £ q in lead I = L,R. As an 
initial condition, we employ the grand-canonical ensemble of the electrons in the leads at 
inverse temperature j8 = \jk&T and with electro-chemical potentials /1l/r- Therefore, the lead 
electrons are described by the equilibrium Fermi function fe(£ q ) = {1 +exp[— fi(£ q — fli))}" . 
For the initial density matrix, we then have 

Pleads,eq « exp [-j8 (H ieads - ^l L N L - ^IrNr)] , (4) 

where N( = Y,q c \ q c iq denotes the electron number in the left and right lead, respectively. 
From this follows that all expectation values of the lead operators can be traced back to the 
expression (c^, ql Cf q ) = S££fd qq ifg(e q ). The two dots couple via the tunnelling matrix element 
Vg q to the state \£q) in the respective lead. The Hamiltonian describing this interaction has the 
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form 

^dot-lead = Y,( V ^ C Lq Cl + V ^ C Rq C2 ^ + HC ' ^ 
1 

It will turn out that the influence of the tunnelling matrix elements is completely characterised 
by the spectral density T((e) = 1%Y*q \Vg q | 2 5(e — e q ) which becomes a continuous function 
of £ if the lead modes are dense. If all relevant lead states are located in the centre of the 
conduction band, the energy-dependence of the spectral densities is not relevant and can be 
replaced by a constant, r L / R (e) = T L / R . This defines the so-called wide-band limit. 

3.2. Master equation approach 

The computation of stationary currents can be achieved by deriving a master equation for the 
dynamics of the dot electrons. Thereby, the central idea is to consider the contact Hamiltonian 
(5) as a pertubation. From the Liouville-von Neumann equation ihp = [H,p] for the total 
density operator p one obtains by standard techniques [18] the approximate equation of 
motion 

i 1 f°° ~ 

Pit) = -r[fldots(0 + Meads, P (03 ~ 71 / dT [^dot-lead, [^dot-lead ( ~ t), p (*)]]• (6) 

n h Jo 

The tilde denotes operators in the interaction picture with respect to the central system and 
the lead Hamiltonian, X(t) = UQ(t)XUo(t), where Uq is the propagator without the coupling. 
The stationary current defined as the net (incoming minus outgoing) electrical current through 
contact £ is given by minus the time-derivative of the electron number in that lead multiplied 
by the electron charge —e, = e(d/dt){Ne). From the master equation (6) follows 

e f°° ~ 

I L (t) = etr[ P (t)N L ] = --2j Q dT<[// dot _ lead (-T), [Mot-lead,^]] )• (7) 

In the following, we specify the master equation (6) and the current formula (7) for 
studying two limiting cases: The first limit U = describes non-interacting electrons. The 
second limit refers to strong Coulomb repulsion such that U is much larger than any other 
energy scale of the problem. Then, only the states with at most one excess electron on the 
wire are relevant. 



3.2.1. Non-interacting electrons In general, the relation between the states |0 a ) and the 
many-particle Hamiltonian (1) is established via the Slater determinant. Alternatively, one 
can resort to Green's functions. In the present case, knowledge of the Green's function at 
time t = is already sufficient. Apart from a prefactor, it is given by the expectation value 
P a p — (cjjC a ) for which one obtains from equation (7) for the stationary current the expression 

eYj 
h 



\(ni\$a)\ 2 fe(£a) 



(8) 



p 

where the index refers to U = 0. It can be shown that the current is independent of the index 
£, i.e. independent of the contact at which it is evaluated. This reflects for a two-probe setting 
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the validity of the continuity equation. For the steady state expectation values P a p, we obtain 
from the master equation (6) the condition 

i(e a -£p)P a p = £ -j{(<l>a\ne)(ni\(l>p)[Mea)+fe(£p)] 



=L,R 



a' p> 
In a non-equilibrium situation, the solution of this set of equations generally possesses non- 
vanishing off-diagonal elements, which in some cases turn out to be crucial. 



3.2.2. Strong Coulomb repulsion In the limit of strong Coulomb repulsion, U is assumed to 
be so large that at most one excess electron resides on the system. Thus, the available Hilbert 
space is restricted to the states {|0),Ca|0)} a= i 5 2, such that the density operator can be written 
as 

p = |0)p o(0|+£(c t a |0)p a o(0| + |0)poa(0|c a )+£4|0)p a/3 (0|c /3 . (10) 

a a/3 
while the current expectation value (7) becomes 

Zoo = er £ £ [£(0j3k)(«£|0a)iKea)paj3 ~ I (0aM \ 2 fe(£a)pOo] , (11) 
a p 

where / = 1 — /. The decomposition of the master equation (6) into the single-particle states 
Ca\0) provides for the stationary state the set of equations 

i(e a -£p)p a p = £ -Y\(<l>a\ne)(ni\<l)p)(fe(£ a )+fe(£p))poo 
£=h,R Z 1 

- £(0a W) (ni\(l>ai)M e a')Pa>p ~ £(0j3' K) ( n Mp )M e p>)Pap> }■ ( 12 ) 
a' p> 

In order to fully determine the density operator, we need in addition an expression for poo 
which can also be derived from the master equation. A more convenient alternative is provided 
by the normalisation condition trp = poo + 'LaPaa = 1- For the sake of completeness, we 
remark that from the master equation (6) follows p a o = poa = in the stationary state. 



3.3. Comparison with experimental data 

Before establishing a quantitative relation between our model and the experimental results, 
we discuss the transport properties of the double-dot model qualitatively. Thereby we reveal 
that both Coulomb repulsion and an orbital degree of freedom play a role for the behaviour 
for fixed, not too small source-drain voltage while the gate voltage is changed. For very 
large negative values of Vgate, both eigenstates of the double dot lie well above the chemical 
potential of both leads and, thus outside the voltage window. This means that lead states 
being in resonance with the dot states remain unoccupied such that electron transport can 
only occur via co-tunnelling processes. Thus, the current will be rather small. When Vg ate 
becomes larger such that the lower dot level lies within the voltage window [see figure 6(a)], 
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Figure 6. Sketch of the tranport through eigenenergy levels for different gate voltages. The 
arrows indicate the possible tunnel evenyts for electrons into and out of of the system. In panel 
(a), VGate is so large that only one level lies within the voltage window, while the other one 
lies well above and is never occupied. Consequently, transport is interaction independent. If 
one level lies below both chemical potentials (b), it will occupied in the steady state and, thus, 
strong Coulomb repulsion inhibits the electron transport. 

resonant transport becomes possible yielding a noticeable current. Increasing the voltage 
further such that also the second level enters the voltage window, opens a second path for 
non-interacting electrons through the dots. In the case of strong Coulomb repulsion, however, 
double occupation of the dot is impossible and, thus, the second orbital cannot fully contribute 
to the transport. Accordingly, the increase of the current is smaller. 

The most significant difference between the two cases is found when only the upper level 
lies within the voltage window, while the lower level is below both chemical potentials, as 
sketched in figure 6(b). Then the stationary state is characterised by an occupied lower level. 
Whether or not a further electron can enter and cause a non- vanishing current now depends on 
the strength of the Coulomb repulsion — for strong repulsion, transport is Coulomb blocked. 
Consequently, for the two limits under investigation, we obtain a current only in the one of 
non-interacting electrons. This is visible as an even qualitative difference in the Coulomb 
diamond structure of figure 3: The scenario for non-interacting electrons complies with 
particle-hole symmetry. This has the consequence that the corresponding Coulomb diamond 
[figure 3(c)] is invariant under changing the sign of both the source-drain voltage and the gate 
voltage. For strong Coulomb repulsion, by contrast, the symmetry concerning the sign of Vg a te 
is no longer present; see figure 3(b). In particular for Vsd ~ 0, the experimental data exhibit 
only one spot with high conductance, which is in clear contrast to the theoretical result for 
non-interacting case shown in figure 3(c). 
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The limits of strong and zero interaction have in common that when the two energy levels 
enter the voltage window from above, the current increases in two steps whose separation is 
determined by the energy splitting 8E of the two levels, i.e. they are separated by the voltage 
8E/e. This enables one to determine the excitation energy of an orbital degree of freedom 
by conductance measurements. Moreover, the source-drain voltage relates to the steepness of 
the triangle. Figure 4 shows again the measurement of figure 3(a) but now with the idealised 
diamond structure marked by dashed lines. 

For a more quantitative treatment, we compare the structure of the experimental result 
in figure 3(a) with the theoretical Coulomb diamond in figure 3(b). This allows one to 
read off the tunnel coupling A = 0.8 meV. For the slightly asymmetric lead-dot couplings 
Tl = 0.2A and I\ = 0.25A, we obtain at the plateaus for the current the values 2.2 nA and 
2.4 nA, respectively, which is of the order of the measured values at the edges of the Coulomb 
diamonds (2 nA — 3 nA). 

The quantitaive agreement between the experiment and the theoretical result for U = °° 
suggests that electrons in the relevant localised states of our sample strongly repell each 
other. This rises the question in which part of the sample (see figure 1) the localised 
states are formed. By the chemical wet etching, the sawtooth pattern has been created and 
four narrow constrictions intersect the long wire into three separate regions. In figure 2, 
features of a one-dimensional system have been detected for less negative side gate voltages 
as discussed in an earlier section. This means that one of the narrower constrictions must 
govern the conduction process in the open channel regime, because it is unlikely, that 
all constrictions represent identical tunnel barriers. However, we cannot identify which 
constriction dominates. Moreover, in samples like the one used in this work, randomly 
distributed charged impurities from the doping process are present. They can strongly 
influence the potential profile depending on their position during the cooling down process 
[19]. 

Unfortunately, without further investigation, we are not able to determine in which part 
of the sample the relevant levels are localised. With the data from the theoretical model for 
strong inter-dot coupling, we nevertheless can infer, that both dots must be rather close. For 
the curved shape of the etched potential into account, a unintentionally emerged dot is very 
likely as well. 

4. Conclusion 

We have studied Coulomb oscillations on lateral fabricated quantum dots near the pinch- 
off. In order to gain a theoretical understanding, we investigated a two-site model which 
implies the consideration of one orbital excitation. A comparison of the measured Coulomb 
diamond with theoretical predictions indicate that both Coulomb repulsion and orbital degrees 
of freedom play a significant role for the transport. The importance of Coulomb repulsion is 
emphasised by the fact that the corresponding model with non-interacting electrons makes 
qualitatively wrong predictions. Moreover, the theoretical results allow us to gauge the gate 
voltage and to determine the energy splitting associated with the orbital excitation. As a 
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drawback, the measurement does yield any conclusion about the nature and the location 
of the two relevant states. Studying a sample in which more orbital degrees of freedom 
play a dominant role might provide additional information. Such experiments should be 
accompannied by theoretical studies for finite Coulomb repulsion strength. 
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